ON THE SPECTRUM OF JACOBI OPERATORS WITH 
QUASI-PERIODIC ALGEBRO-GEOMETRIC COEFFICIENTS 



VLADIMIR BATCHENKO AND FRITZ GESZTESY 



Abstract. We characterize the spectrum of one-dimensional Jacobi operators 
H = aS + + a~ S~ + b in i 2 (Z) with quasi-periodic complex- valued algebro- 
geometric coefficients (which satisfy one (and hence infinitely many) equa- 
tion^) of the stationary Toda hierarchy) associated with nonsingular hyperel- 
liptic curves. The spectrum of H coincides with the conditional stability set 
of H and can explicitly be described in terms of the mean value of the Green's 
function of H. 

As a result, the spectrum of H consists of finitely many simple analytic 
arcs in the complex plane. Crossings as well as confluences of spectral arcs are 
possible and discussed as well. 



1. Introduction 

It is well-known from the work of Date and Tanaka [12] , [13] , Dubrovin, Matveev, 
and Novikov [17], Flaschka [23], McKean [40], McKean and van Moerbeke [41], [42], 
Mumford [43], Novikov, Manakov, Pitaevskii, and Zakharov [47], Teschl [53, Chs. 
9,13], Toda [54, Ch. 4], [55, Chs. 26-30], van Moerbeke [56], van Moerbeke and 
Mumford [57], that the self-adjoint Jacobi operator 

H = aS + + a~S~ +b, dom(H) = £ 2 (Z), (1.1) 

in £ 2 (Z) with real- valued periodic, or more generally, algebro-geometric quasi- 
periodic and real-valued coefficients a and b (i.e., coefficients that satisfy one (and 
hence infinitely many) equation(s) of the stationary Toda hierarchy), leads to a 
finite-gap, or perhaps more appropriately, to a finite-band spectrum o~(H) of the 
form 

P+i 

a{H) = (J [E 2 

in — 

m— 1 

Compared to the real-valued case, the corresponding spectral properties of Jacobi 
operators with periodic and complex-valued coefficients a and b, to the best of our 
knowledge, have been studied rather sparingly in the literature. We are only aware 
of two papers by Naiman [45] , [46] , in which it is shown that the spectrum consists 
of a set of piecewise analytic arcs which may have common endpoints. (This is 
in anology to the case of (non-self-adjoint) one-dimensional periodic Schrodinger 
operators, cf. [48], [51].) It seems plausible that the latter case is connected with 
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(complex- valued) stationary solutions of equations of the Toda hierarchy. In partic- 
ular, the next scenario in line, the determination of the spectrum of H in the case of 
quasi-periodic and complex-valued solutions of the stationary Toda hierarchy appar- 
ently has never been clarified. The latter spectral problems have been open since 
the late seventies and it is the purpose of this paper to provide a comprehensive so- 
lution of them. For theta function representations of a and b in the complex- valued 
algebro-geometric case (without addressing spectral theoretic questions) we refer, 
for instance, to Aptekarev [1], Dubrovin, Krichever, and Novikov [16], Krichever 
[33]-[36] (cf. also the appendix written by Krichever in [15]). 
To describe our results, a bit of preparation is needed. Let 

G{z, n, ri) = {H- z)~\n, n'), z e <C\o{H), n, n' e Z, (1.3) 

be the Green's function of H (here cr(H) denotes the spectrum of H) and denote 
by g(z,n) the corresponding diagonal Green's function of H defined by 

g(z, n) = G(z, n, n) = , , (1.4) 

n>2p+2\z) i 

2p+l 

R2 P+ 2{z) = J] ( z - E ™)> {^>m=o C C, (1.5) 
E m ^E m > for m ^ m', to, m! = 0, 1, . . . , 2p + 1. (1.6) 

For any quasi-periodic (in fact, Bohr (uniformly) almost periodic) sequence / = 
{/(^)}feez the mean value (/} of / is defined by 
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N 



(f) = lim n ^ i V* f( k )- (1-7) 

k=-N 



Moreover, we introduce the set S by 
S = | A G C 



where y 2 = R.2p+2{z), and the polynomial G p +\(z, n) of degree p + 1 in z will be 
defined in (2.21). Here we observe that G p +i is given in terms of the off-diagonal 
Green's function G(z,n+ l,n) by 

Gp+l(z ' n) -l-2a(n)G(z,n+l,n) (1.9) 



R2 P +2{z) 1 /'- 
(cf. also (3.13)). In addition, we note that 

XUy'* (lio) 

for some constants {Aj}^ =1 C C. 

Finally, we denote by a p (T), a T (T), a c (T), a e (T), and cr ap (T), the point spectrum 
(i.e., the set of eigenvalues), the residual spectrum, the continuous spectrum, the 
essential spectrum (cf. (4.16)), and the approximate point spectrum of a densely 
defined closed operator T in a complex Hilbcrt space, respectively. 

Our principal new results, to be proved in Section 4, then read as follows: 
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Theorem 1.1. Assume that a and b are quasi-periodic (complex-valued) solu- 
tions of the pth stationary Toda equation associated with the hyperelliptic curve 
y 2 = i?2p+2(z) subject to (1.5) and (1.6). Then the following assertions hold: 

(i) The point spectrum and residual spectrum of H are empty and hence the spec- 
trum of H is purely continuous, 

<t p (H) = * i (H) = Q, (1.11) 
a(H) = a c (H) = a e (H) = a ap (H). (1.12) 

(ii) The spectrum of H coincides with £ and equals the conditional stability set of 



a(H) = Ae. ... 

\ V\ \G p+1 (\-)+ V/ 

= {A G C | there exists at least one bounded solution 

^ V e f°(Z) of Hip = (1.14) 

(Hi) a(H) C C is bounded, 

o-(H) C {z G C | Re(z) G [Mi, M 2 ], hn(z) G [M 3 , M 4 ]}, (1.15) 



where 



Mi = -2sup[|Re(a(n))|] + inf [Re(6(n))], 



M 2 = 2sup[|Re(a(n))|] +sup[Re(6(n))], 



/it: 



M 3 = -2sup[|Im(a(n))|] + inf pm(6(n))], 
M 4 = 2sup[|Im(a(n))|] + sup[Im(6(n))]. 

n£Z neZ 



(1.16) 



(it)) cr(i?) consists of finitely many simple analytic arcs. These analytic arcs may 
only end at the points Ai, . . . , X p , Eq, . . . , i?2p+i • 

(v) Each E m , m = 0, . . . , 2p+l, is met by at least one of these arcs. More precisely, 
a particular E mo is hit by precisely 2Nq + 1 analytic arcs, where Nq G {0, . . . ,p} 
denotes the number of Xj that coincide with E mo . Adjacent arcs meet at an angle 
2n/(2No + 1) at E mo . (Thus, generically, N — and precisely one arc hits E mo .) 

(vi) Crossings of spectral arcs are permitted and take place precisely when 

Re(( ln(' Gp+l( ^°' ,) ~ V \))=0 



,G P+1 (\ JO ,-) + yJ / J ^ (1.17) 
for some j G {1, . . . ,p} with X jo £ {E m } 2 ^Xl- 

In this case 2M + 2 analytic arcs are converging toward \j , where M G {!,... ,p\ 
denotes the number of Xj that coincide with Xj . Adjacent arcs meet at an angle 
7r/(M + 1) at Xj . (Thus, if crossings occur, generically, M = 1 and two arcs 
cross at a right angle.) 

(vii) The resolvent set C\a(H) of H is path- connected. 

Naturally, Theorem 1.1 applies to the special case where a and b are periodic 
complex-valued solutions of the pth stationary Toda equation associated with a 
nonsingular hyperelliptic curve. Even in this special case, Theorem 1.1 yields new 
facts which go beyond the previous results by Naiman [45] , [46] . 
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For analogous results in the context of one-dimensional Schrodingcr operators 
with quasi-periodic algebro-geometric KdV potentials we refer to [2] . 

Theorem 1.1 focuses on stationary quasi-periodic solutions of the Toda hierarchy 
for the following reasons. First of all, the class of algebro-geometric solutions of 
the (time-dependent) Toda hierarchy is defined as the class of all solutions of some 
(and hence infinitely many) equations of the stationary Toda hierarchy. Secondly, 
time-dependent algebro-geometric solutions of a particular equation of the (time- 
dependent) Toda hierarchy just represent isospectral deformations (the deformation 
parameter being the time variable) of fixed stationary algebro-geometric Toda so- 
lutions (the latter can be viewed as the initial condition at a fixed time to)- In the 
present case of quasi-periodic algebro-geometric solutions of the pth Toda equation, 
the isospectral manifold of such given solutions is a complex p-dimensional torus, 
and time-dependent solutions trace out a path in that isospectral torus (cf. the 
discussions in [24] and [25]). 

Finally, we give a brief discussion of the contents of each section. In Section 2 
we provide the necessary background material including a quick construction of the 
Toda hierarchy of nonlinear evolution equations and its Lax pairs using a polynomial 
recursion formalism. We also discuss the hyperelliptic Riemann surface underlying 
the stationary Toda hierarchy, the corresponding Baker-Akhiezer function, and 
the necessary ingredients to describe the analog of the Its-Matveev formula for 
stationary Toda solutions. Section 3 focuses on the Green's function of the Jacobi 
operator H, a key ingredient in our characterization of the spectrum cr(H) of H 
in Section 4 (cf. (1.13)). Our principal Section 4 is then devoted to a proof of 
Theorem 1.1. Appendix A provides the necessary summary of tools needed from 
elementary algebraic geometry (most notably the theory of compact (hyperelliptic) 
Riemann surfaces) and sets the stage for some of the notation used in Sections 2-4. 
Appendix B provides additional insight into one ingredient of the theta function 
representation of the coefficients a and b. 

2. The Toda Hierarchy, Hyperelliptic Curves, and Theta Function 
Representations of the Coefficients a and b 

In this section we briefly review the recursive construction of the Toda hierarchy 
and associated Lax pairs following [7] and [25]. Moreover, we discuss the class of 
algebro-geometric solutions of the Toda hierarchy corresponding to the underlying 
hyperelliptic curve and recall the analog of the Its-Matveev formula for such so- 
lutions. The material in this preparatory section is known and detailed accounts 
with proofs can be found, for instance, in [7]. For the notation employed in connec- 
tion with elementary concepts in algebraic geometry (more precisely, the theory of 
compact Riemann surfaces), we refer to Appendix A. 

Throughout this section we assume that 

a, b e r°(Z), a(n) ^ for all n G Z, (2.1) 
and consider the second-order Jacobi difference expression 

L = aS+ + a-S- +b, (2.2) 
where denote the shift operators 



(S ± f)(n) = /±(n) = /(nil), n e Z, / e i°°(Z). (2.3) 
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To construct the stationary Toda hierarchy we need a second difference expres- 
sion of order 2p + 2, p G No, defined recursively in the following. We take the 
quickest route to the construction of P 2p + 2 , and hence to the Toda hierarchy, by 
starting from the recursion relations (2.4)-(2.6) below. 

Define {fj}jen and {gj}jeN recursively by 

/o = l, go = -c u (2.4) 

2f j+1 + gj + gj - 2b fj = 0, jeN , (2.5) 

g j+ i-g; +1 + 2(a 2 ff -(a-) 2 ff)-b( 9j -gr) = 0, jeN„. (2.6) 
Explicitly, one finds 

/o = l, 
fi =b + a, 

f 2 = a 2 + (a~) 2 + b 2 + Cl b + c 2 , etc., (2.7) 

.9o = -ci, 

gx = -2a 2 - c 2 , 

g 2 = -2a 2 (b + + b)+ Cl (-2a 2 )-c 3 , etc. 

Here {cj}j e ^ denote undetermined summation constants which naturally arise when 
solving (2.4)-(2.6). 

Subsequently it will be convenient to also introduce the corresponding homoge- 
neous coefficients fj and gj, defined by vanishing of the constants c k , k G N, 

/0 = 1. fi = fi\ Ch =0,k=l,...J' -?' GN ' 

9J = 9i\ Ck=0tk=l _ j+v ieN . (2.8) 



Hence. 



introducing 



/j = ^2 c i-kfk, gj = ^2 c o-kh - Cj+1, j G N , (2.9) 

fe=0 fe=l 



c = 1. (2.10) 
Next we define difference expressions P 2p +2 of order 2p + 2 by 

v 

-LP+ 1 + J2 (ft + 2o/ J -5+)LP-j + /p+i, p G N . (2.11) 

3=0 

Using the recursion relations (2.4)-(2.6), the commutator of P 2p + 2 and L can be 
explicitly computed and one obtains 

[P 2p+2 , L] = - a(g+ + g p + f+ +1 + f p+1 - 2b+f+)S+ 

+ 2( - b(g p + f p+1 ) + a 2 f+ - {a-ff- + b 2 f p ) 

- ^ (g P + g~ + f P+ i + f~ +1 - 2bf p )S~, p G N . (2.12) 

In particular, {L,P 2p+2 ) represents the celebrated Lax pair of the Toda hierarchy. 
Varying p G No, the stationary Toda hierarchy is then defined in terms of the 
vanishing of the commutator of P 2p + 2 and L in (2.12) by, 

[P 2p+2 , L] = s-Tl p (a, b) = 0, p G N . (2.13) 
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Thus one finds 

9p + 9p + f P +i + f~+i - 26/ P - 0, (2.14) 

-b(g P + U+i) + a 2 f+ - (a-ffg + b 2 f P = 0. (2.15) 
Using (2.5) with j = p one concludes that (2.14) reduces to 

f P +i = / P +i, (2-16) 

that is, f p+ i is a lattice constant. Similarly, one infers by subtracting b times (2.14) 
from twice (2.15) and using (2.6) with j — p, that g p+ i is a lattice constant as well, 
that is, 

9 P +i=9p+i- (2-17) 

Equations (2.16) and (2.17) constitute the pth stationary equation in the Toda 
hierarchy, which will be denoted by 

s-Tb>, b) = ( ~ fp _ +1 ) = 0, pe No. (2.18) 
\9 P +i -9p+\J 

Explicitly, 



s-Tl (a,6) = 
s-Tli(a,6) = 




(2.19) 



represent the first few equations of the stationary Toda hierarchy. By definition, 
the set of solutions of (2.13), with p ranging in No and Cfe G C, k € N, represents 
the class of algebro-geometric Toda solutions. 

In the following we will frequently assume that a and b satisfy the pth stationary 
Toda equation. By this we mean they satisfy one of the pth stationary Toda 
equations after a particular choice of integration constants Ck <E C, k — l,...,p, 
peN, has been made. 

Next, we introduce polynomials F p (z,n) and G p +i(z,n) of degree p and p + 1 
with respect to the spectral parameter z € C by 

p 

Fp{z,n) = Y J Z i fp-j{n), (2.20) 

3=0 

V 

G p+1 (z,n) = -z p+1 + J2 zj 9 P -j(n) + f P +i(n). (2.21) 
j=o 

Explicitly, one obtains 
F = l, 

F 1 =z + b + c 1 , 

F 2 = z 2 + bz + a 2 + (a-) 2 +b 2 + Cl (z + b) + c 2 , etc., (2.22) 
d = -z + b, 

G 2 = -z 2 + (gT) 2 - a 2 + b 2 + d(-z + b), etc. 
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Next, we study the restriction of the difference expression P2 P +2 to the two- 
dimensional kernel (i.e., the formal null space in an algebraic sense as opposed to 
the functional analytic one) of (L — z). More precisely, let 

kcr(L -z) = {y;:Z^CU {00} | (L - z)i\) = 0}. (2.23) 
Then (2.11) implies 

P2 P+ 2 \ kel{L - z) = (2aF p (z)S + + G p+1 (z))\ kci(L _ zy (2.24) 
Therefore, the Lax relation (2.12) becomes 

= [P 2p+2 , L] | kor( L-,)= (a(2(z - b+)F+ - 2{z - b)F p + G~ +1 - G+ +1 )S+ 

+ (2(a-) 2 F- - 2a 2 F+ + (z - b)(G~ +1 G p+1 ))) j^^, (2.25) 
or, equivalently, 

2(z - b+)F+ - 2{z - b)F p + G+ +1 - G~ +1 = 0, (2.26) 

2a 2 F+ - 2(a-) 2 F- + (z b)(G p+1 G~ +1 ) = 0. (2.27) 
Upon summing (2.26) one infers 

2(z - b+)F+ + G+ +1 + G p+1 = 0, peN , (2.28) 
and inserting (2.26) into (2.27) then implies 

(z - b) 2 F p + {z- b)G p+1 + a 2 F+ - (a-) 2 F p - = 0, P e N . (2.29) 
Combining equations (2.27) and (2.28) one concludes that the quantity 

R2 P+2 (z) = G p+1 (z,n) 2 - Aa{n) 2 F p {z,n)F+{z in ) (2.30) 
is a lattice constant, and hence depends on z only. Thus, one can write 

2p+l 

R2 P+ 2(z) = 1] ( z - E ™)> { E m?Ztl C C. (2.31) 

One can show that equation (2.30) leads to an explicit determination of the inte- 
gration constants c\ , . . . , c p in 

s-Tl p (a,6) = (2.32) 

in terms of the zeros E , . . . , E 2p +i of the associated polynomial R 2p +2 in (2.31). 
In fact, one can prove that 

c k = c k (E), k = l,...,p, (2.33) 

where 



-J2p + 1 



k{ ~ ] ~ jo j^, 2 2fe (jo!) 2 • • • (J 2p+ i!) 2 (2. 70 -!)••• (2h P+1 if 

joH h72p+i = k 



(2. 7 o)!---(2j 2 p+i)!_ 

• (2j s 

k = l,...,p (2.34) 

are symmetric functions of E_ = (Eq, . . . , -E2p+i). 

Remark 2.1. Since by (2.5), (2.6), (2.20) and (2.21), a enters quadratically in F p 
and G p+ i, the Toda hierarchy (2.18) is invariant under the substitution 

a(n) a e (n) = e(n)a(n), e(n) e {+1, -1}, n e Z. (2.35) 
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We emphasize that the result (2.24) is valid independently of whether or not 
P2 P +2 and L commute. However, the fact that the two difference expressions P2 P +2 
and L commute implies the existence of an algebraic relationship between them. 
This gives rise to the Burchnall-Chaundy polynomial for the Toda hierarchy first 
discussed in the discrete context by Naiman [45] , [46] . 

Theorem 2.2. Fixp G No and assume that P2 P +2 and L commute, [P2 P +2, L] = 0, 
or equivalently, suppose s-Tl p (a, b) = 0. Then L and ?2p+2 satisfy an algebraic 
relationship of the type (cf. (2.31)) 

Fp(L, P2P+2) = Pip+2 - R2p+2{L) = 0, 

^ (2.36) 
R2p+2{z)= Yl(z-E m ), ze<C. 

The expression T P {L, P2P+2) is called the Burchnall-Chaundy polynomial of the 
Lax pair (L,P 2 p+2)- Equation (2.36) naturally leads to the hyperelliptic curve fC p 
of (arithmetic) genus p€No, where 

K p : T p (z, y)=y 2 - R 2 p+2{z) = 0, 

2p+l / 2 37) 

R2 P+ 2(z) =\{{z E m ), {E m fZtl C C. 

m=0 

The curve JC P is compactified by joining two points Poo ± , Poo + ^ -Poo_ , at infinity. 
For notational simplicity, the resulting curve is still denoted by IC p . Points P on 
/C p \P 00± are represented as pairs P — (z, y), where y(-) is the meromorphic function 
on Kp satisfying !F p (z,y) = 0. The complex structure on JC p is then defined in the 
usual way, see Appendix A for the case of nonsingular curves. Hence, IC P becomes 
a two-sheeted hyperelliptic Riemann surface of (arithmetic) genus p € No (possibly 
with a singular affine part) in a standard manner. 

We also emphasize that by fixing the curve JC P (i.e., by fixing Eq, . . . ,i?2p+i), 
the integration constants c\ , . . . , c p in the corresponding stationary s-Tl p equation 
are uniquely determined as is clear from (2.33) and (2.34), which establish the 
integration constants Ck as symmetric functions of Eq, . . . , E2 P +\. 

For notational simplicity we will usually tacitly assume that p E N. The trivial 
case p = 0, which leads to a(n) 2 = {E 1 - E ) 2 /16, b(n) = (E + E^/2, n e N, is of 
no interest to us in this paper. 

In the following, the zeros 1 of the polynomial F p (-,n) (cf. (2.20)) will play a 
special role. We denote them by {^j{n)} p - =1 and write 

p 

F p (z,n) = Y[[z - ixjin)]. (2.38) 

i=i 

The next step is crucial; it permits us to "lift" the zeros fMj of F p from C to the 
curve IC p . From (2.30) and (2.38) one infers 

R2 P +2(z) - G p+ i{zf =0, z e {Mj.Mfe}j,fe=i,...,p- ( 2 -39) 
We now introduce {Aj( n )}j=i,...,p C JC P by 

p,j(n) = (fj,j(n),-G p+1 (fj,j(n),n)), j = 1, ...,p, n G Z. (2.40) 

1 If a, b g £°°(Z), these zeros are the Dirichlet eigenvalues of a bounded operator on £ 2 (Z) 
associated with the difference expression L and a Dirichlet boundary condition at n 6 Z. 
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Next, we recall equation (2.30) and define the fundamental meromorphic function 
<f>(-, n) on JCp by 

flp,n) = ^%fag (2-41) 
^ y ' ; 2a(n)F p (z,n) y ' 



-2a(n)F p (z,n+l) 



(2.42) 



V + G p+1 (z,n) 
P = (z, y) e Kp, neZ 

with divisor {<j>{-, n)) of </>(•, n) given by 

(0(-,fl)) - Pp oc+ A(n+l) - ^Poo_ AW' ( 2 ' 43 ) 

using (2.38) and (2.40). Here we abbreviated 

A = {Ai,---,Ap}eSymf(/C p ) (2.44) 

(cf. the notation introduced in Appendix A) . The stationary Baker- Akhiezer func- 
tion ip(-, n, no) on KL P is then defined in terms of (/)(■, n) by 

nrino^^) for «>«o+i, 

ip(P, n,n ) = < 1 for n = n , (2.45) 
.nr^^m)" 1 for n<n -l 

with divisor (^(-, n, no)) of 4>(P, n, no) given by 

(V(-,n,no)) = %„) - %„ o) + (n - no)(X> Poo+ - X> Poo _). (2.46) 

Basic properties of (f> and -0 are summarized in the following result. We denote by 
ff)(n) = a(fg + — f + g) the Wronskian of two complex- valued sequences / and 
g, and denote P* = (z, —y) for P = (z, y) e /C p . 

Lemma 2.3. Assume (2.1) and suppose a,b satisfy the pth stationary Toda equa- 
tion (2.18). Moreover, let P = (z,y) e /C P \{Poo±} and (n, no) G Z 2 . TTien 
satisfies the Riccati-type equation 

acp(P) +a-<p-(P)- 1 = z-b, (2.47) 

as well as 

*(P)MP") = (248) 

^^^Sgir < 2 - 49 » 

Moreover, ip satisfies 

(L — z(P))tl>(P) = 0, (P 2p+2 -y(P))^(P) = 0, (2.51) 

^(P, n, n )^(P*, n, n ) - t^t^t, (2.52) 

F p (z,n ) 



a(n) [ip(P, n, n )ip(P*,n + 1, n ) + i}j{P*,n, n )ip(P, n + 1, n Q )] 
= G p+ i(z,n) 
P p (z,n ) 



(2.53) 
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WMP, ■, n ), ^(P*, •, n )) = - F V , { p n y (2.54) 

rp{z, no) 

Combining the polynomial recursion approach with (2.38) readily yields trace 
formulas for the Toda invariants, which are polynomial expressions of a and b, in 
terms of the zeros /Uj of F p . 

Lemma 2.4. Assume (2.1) and suppose a,b satisfy the pth stationary Toda equa- 
tion (2.18). Then, 

°(») 2 = ^E^+ 2 (A») f[[^(n) - /ife(n)]- 1 + l[6( 2 )(n) - &(„)*], (2.55) 
j=i fe=i 



fe H = ^E^-E^W, (2.56) 

m— j — 1 

2p+l p 

& (fc) (n)=2E^-E^(") fc > fceR ( 2 ' 5? ) 

m— j=l 

Strictly speaking, (2.55) is only valid in the case where for all n G Z, /ij(n) ^ 
Ufc(n) for j 7^ fc. The case where some of the fij coincide requires a limiting 
argument that will be omitted as (2.55) will play no further role in this paper. The 
details of this limiting procedure can be found in [26] . 

From this point on we assume that the afhne part of JC P is nonsingular, that is, 

E m ^ E m , for m^m', m, to' = 0, 1, . . . , 2p + 1. (2.58) 

Since nonspecial divisors play a fundamental role in this context we also recall 
the following fact. 

Lemma 2.5. Assume (2.1) and suppose a,b satisfy the pth stationary Toda equa- 
tion (2.18). In addition, assume that the affine part of K p is nonsingular. Let V^, 
/i = {/ii, ... ,/tp} G Sym p (£p), be the Dirichlet divisor of degree p associated with 
a, b defined according to (2.40), that is, 

(tj(n) = (itj(n),-G p+1 (fij(ri),n)), j = 1, . . . ,p, n G Z. (2.59) 

Then T^jx( n ) * s nonspecial for all n G Z. Moreover, there exists a constant C M > 
such that 

Mn)| < C„, j = l,...,p, nGZ. (2.60) 

We remark that if a, b G £°°(Z) satisfy the pth stationary Toda equation (2.18), 
then automatically a(n) ^ for all n G Z (cf. [26]). 

We continue with the theta function representation for ip, a, and b. For general 
background information and the notation employed we refer to Appendix A. 

Let 9 denote the Riemann theta function associated with K p (whose affine part is 
assumed to be nonsingular) and a fixed homology basis {dj, bj}^ =1 . Next, choosing 
a base point Qq G B(K. p ) in the set of branch points of 1C P , we recall that the Abel 
maps Aq q and o_q () are defined by (A. 40) and (A. 43), and the Riemann vector Sq o 
is given by (A. 55). Then Abel's theorem (cf. (A. 53)) and (2.46) yields 

^Qo (%«))= S Qo (%«„)) - A Poo (Poo + )(n - n ) 

(2.61) 

= «q,Paw) - 2A Qo (P oc+ )(n - n ). 
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(3) 

Next, let uj-p p denote the normalized differential of the third kind defined 

'C» _J_ ) -* OO _ 

by 



x,.V „ ~\\ -V"/- :; : (s ' • 0(1)H as /' ■ 1\ .. (2.62) 



/3) 



C = V*, 

where the constants Xj e C, j = 1, . . . ,p, are determined by employing the normal- 
ization 

' "pLp^ =0, J = 1,...,P. (2-63) 



/ 



One then infers 



(3) _ J-W^ i J 3 ) 



P^P 



± ln(C) + e " J (Q ) + 0(C) as P -» P^ (2.64) 

Qo 

for some constant e 3 ^(Qo) € C. The vector of fe-pcriods of ujp^ p ^ /(27n) is 
denoted by 

ui 3) = (u!$,...,U$), rZ ^-[ r ;l.r^ ■ 3 = 1,..., P. (2-65) 

Since Qo is a branch point, Qo G &(fcp), one concludes by (A. 42) that 

= A Poo _(P 00+ ) = 2A Qo (P 00+ ). (2.66) 

In the following it will be convenient to introduce the abbreviation 

z(P, Q) = S Qo - A Qo (P) + a Qo (Vq), (2.67) 
P G JCp, Q = {Qi, ...,Q P }e Sym p (/C p ). 

We note that z(-,Q) is independent of the choice of base point Qo- 

The zeros and the poles of ip as recorded in (2.46) suggest consideration of the 
following expression involving ^-functions (cf. (A. 29)) 

*fa(P.gM)) K/<,/^--)- (2 ' 68) 

Here we agree to use the same path of integration from Q to P on JC P in the 
Abel map Aq„ (P) m z{P, A( n )) and in the integral of lo p p in the exponent of 
(2.68). With this convention the expression (2.68) is well-defined on JC P (cf. Remark 
A. 4, however) and one concludes 

m n,n ) = C{n, noj ^g^ exp [(n - n ) ^ o, Poo+iPoe _ j . (2.69) 

To determine C(n, n ) one can use (2.52) for P = P^ and P* = Poo_ . Hence, 
r( , 2 _ -9(^00+. A(no)))gU(Poc + ,A(no - 1))) 

° ln ' noj ~ ^(P^.AWJjtfU^co+.A^-i))) ■ ( j 

Thus, one obtains the following well-known result. 
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Theorem 2.6. Assume (2.1) and suppose a, b satisfy the pth stationary Toda equa- 
tion (2.18). In addition, assume the affine part of K v to be nonsingular and let 
P G /C p \{Poo ± } and n, n G Z. Then T^jj,( n ) * s nonspecial for n G 7L. Moreover, 2 



fl(z(P,A(n))) / ^ (3) 

VOP, », n ) = C(n, "o) gU(p ^ (nq))) exp ^(n - no) ^ 



where 



C(n, n ) 



%(P>o + ,A(rio)))%(P>o + ,A>o - !))) 



1/2 



%(Poo + ,M«)))%(^oo + ,£(™ - 1))) 

w«i/i i/ie linearizing property of the Abel map, 

Sq (%«)) = (fiQ (^A(«o)) - 2A Qo (P 00+ )(n - n )) (mod L p ). 
The coefficients a and b are given by 

"%(Poo + ,A(n - l)))0(z(P oo+ ,A(ra+ 1)))1 1/2 



a(n) = a 



0(«(Poo+, £(«))) 
K») = -z E E ™ - E A -> + E c » 7TT ln 



n G 



(2.71) 

(2.72) 
(2.73) 

(2.74) 



2p+l 



m=0 



+ z(p oc+ ,A(n))) 



w=0 

(2.75) 



_%+z(P 00+ ,A>-l))) 
n G Z 

where the constant a depends only on JC P and cj(p) is given by (A. 23). 

Combining (2.73) and (2.75), one observes the remarkable linearity of the theta 
function with respect to n in formulas (2.74), (2.75). In fact, one can rewrite (2.75) 
as 



b(n) 



a V- / x 9 , (6(w + A-Bn)\ 



3 = 1 



w=0 



where 



A = S Qo - Aq q (POO+) + U^n + a Qo (%„ o) ), 
C = A + B, 

2p+l 



A 



2p+l p 

Je^-e* 



(2.76) 

(2.77) 

(2.78) 
(2.79) 

(2.80) 



m— 



Hence, the constants A G C and P G C p are uniquely determined by JC p (and its 
homology basis), and the constant A G C p is in one-to-one correspondence with the 
Dirichlet data A( n o) = (Ai( n o), • • ■ , Ap("o)) G Sym p /C p at the point n . 

Remark 2.7. If one assumes a, 6 in (2.74) and (2.75) to be quasi-periodic (cf. (3.6) 

~ ~(3) 

and (3.7)), then there exists a homology basis {aj, bj}*j =1 on K, p such that B_=U 
satisfies the constraint 

l = fl 3) €P. (2.81) 
This is studied in detail in Appendix B. 



2 To avoid multi-valued expressions in formulas such as (2.71), etc., we agree to always choose 
the same path of integration connecting Qo and P and refer to Remark A. 4 for additional tacitly 
assumed conventions. 
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3. The Green's Function of H 

In this section we focus on the properties of the Green's function of H and derive 
a variety of results to be used in our principal Section 4. 
Introducing 



G(P,m,n) = 



1 \ tp(P* ,m,no)ip(P,n,no), m<n, 

W(^{P*, na),il>(P, n )) \i>(P, m, n )^(P*, n, n ), m > n, 

PeWocJ, n,n GZ, (3.1) 

and 

g(P,n) = G(P,n,n) = ^n^P* n,n ) 

equations (2.52) and (2.54) then imply 

9(P,n) = -^^, P=(z,y)elC p \{P 00± }, neZ. (3.3) 

Together with g(P,n) we also introduce its two branches g±(z, n) defined on the 
upper and lower sheets H± of K p (cf. (A. 3), (A. 4), and (A. 14)) 

F (z n) 

9±(z,n) = T p P ; / 1/2 , zeU,neZ (3.4) 

with II = C\C the cut plane introduced in (A. 4). 

For convenience we shall focus on g_ whenever possible and use the simplified 
notation 

g(z,n) = g~(z, n), z£U 7 neZ (3.5) 

from now on. 

Next we briefly review a few properties of quasi-periodic and almost-periodic 
discrete functions. 

We denote by QP(Z) and AP(Z) the sets of quasi-periodic and almost periodic 
sequences on Z, respectively. 

In particular, a sequence / is called quasi-periodic with fundamental periods 
(ill, . . . , Ojv) £ (0, oo)^ if the frequencies 2-ir/fli, . . . , 2ir/£ljy are linearly indepen- 
dent over Q and if there exists a continuous function F e C(R N ), periodic of period 
1 in each of its arguments, 

F(xi,...,x j + l,...,x N )=F(x 1 ,...,x N ), xj e E, j = 1, . . . ,N, (3.6) 

such that 

f(n) =F(fl^ 1 n,...,n^ 1 n), n £ Z. (3.7) 

Any quasi-periodic sequence on Z is almost periodic on Z. Moreover, a sequence 
/ = {f(k)}kei, is almost periodic on Z if and only if there exists a Bohr almost 
periodic function g on R such that f(k) = g{k) for all fc £ Z (see, e.g., [11, p. 47]). 

For any almost periodic sequence / = {.f(k)}kez, the mean value (/) of /, 
defined by 

, n +N 

(f) = lim — — - V f(k), (3.8) 

k=n -N 

exists and is independent of n £ Z. Moreover, we recall the following facts for 
almost periodic sequences that can be deduced from corresponding properties of 
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Bohr almost periodic functions, see, for instance, [5, Ch. I], [6, Sects. 39-92], [11, 
Ch. I], [22, Chs. 1,3,6], [30], [38, Chs. 1,2,6], and [50]. 

Theorem 3.1. Assume f, g G AP(Z) and n n ,n eZ. Then the following assertions 
hold: 

(i) f_et°°(Z). 

(ii) f, cf, ceC, /(■ + n), f(n ■), n G Z, \f\ a , a > are all in AP(Z). 
(in) / + 5 ,/seAP(Z). 

(iv) 1/g G AP(Z) if and only ifl/g G £°°(Z). 

(w) Let G be uniformly continuous on M C R and /(n) G A4 /or a/Z neZ. JTien 

G(.f) G AP(Z). 
(w) Let (/) = 0, ften EL„ , = °(H)- 

| n | — » oc 

(vii) Let F(n) = Y2=n f( k )- Then F e AP ( Z ) */ and onl V l f F e ^°°(^)- 
(«t) J/0 < / G AP(Z)° / ^ 0, then (/) > 0. 

(iar) If l/f G ^°°(Z) and / = |/| exp(^), tfien |/| G AP(Z) and <p is of the type 
tp(n) = cn + tp(n), where c G R and ip G AP(Z) (and real-valued) . 

(x) IfF(n) = cxp(j2l =no /(*;)), then F G AP(Z) if and only if f(n) = i/3 + i>(n), 

where [3 G R, tp G AP(Z), and * G ^°°(Z), w/iere *(n) = ELn ^( fc )- 

For the rest of this paper it will be convenient to introduce the following hypoth- 
esis: 

Hypothesis 3.2. Assume the affine part of /C p to be nonsingular. Moreover, 
suppose that a, b G QP(Z) satisfy the pth stationary Toda equation (2.18) on Z. 

Next, we note the following result. 

Lemma 3.3. Assume Hypothesis 3.2. Then all z-derivatives of F p (z, •) and 
G p +i(z, ■), z G C, and g(z 7 ■), z G II, are quasi-periodic. Moreover, taking limits to 
points on C, the last result extends to either side of cuts in the set C\{E m }^^ (cf. 
(A. 3)) by continuity with respect to z. 

Proof. Since ft and gi are polynomials with respect to a and 6, ft and gt, £ € N, 
are quasi-periodic by Theorem 3.1. The corresponding assertion for F p (z, •) then 
follows from (2.20) and that for g(z, •) follows from (3.4). □ 

In the following we represent G p+ i(z,n) + G p+1 (z,n) as 

p+1 

G p+1 (z,n) + G+ +1 (z,n) = -2 J[[z - Vk(n)], z G C, n G Z, (3.9) 

fe=i 

and note that the roots are bounded, 

IK||oo<<5, fc = l,...,p+l (3.10) 

for some constant C > 0, since the coefficients of G p+ \(z,n) are defined in terms 
of bounded coefficients a and b by(2.6). For future purposes we introduce the set 

n c = n\({zGC||z| <c + i}u 

[z G C I min [Re(E m )] - 1 < Re(z) < max [Re(£ ro )] + 1, 

L m=0,...,2p+l m=0,...,2p+l 

min [Im(£ m )] - 1 < lm(z) < max [Im(£ m )] + l}) , (3.11) 

m=0,...,2p+l m=0,...,2p+l J / 
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where C = max{C M , ||6||oo,C} and C M is the constant in (2.60). Without loss of 
generality, we may also assume that He contains no cuts, that is, 

n c nc = 0. (3.12) 

Next, we derive a fundamental equation for the mean value of the diagonal 
Green's function g(z, •) that will allow us to analyze the spectrum of the Jacobi 
operator H. First, we note that by (2.49), (2.50), (2.54), and (3.1) one obtains 



G(P,n,n+l) _ G p+1 (z,n) - y 



, P = (z,y) e JC p , n £ 



(3.13) 



G(P*,n,n+l) G p+1 (z,n)+y' * v>»/-" <p> 

Differentiating the logarithm of the expression on the right-hand side of (3.13) with 
respect to z and using (2.30), one infers 



2 dz^ n \G p+ i(z,n) + y 



1 d fG p+1 (z,n)-y\ _ -^^G p+1 (z,n)-yG; +1 (z,n) 



, z g n 



c- 



(3.14) 



-Aa(n) 2 F p (z, n)Fp(z, n) 

Here • abbreviates d/dz. We note that the left-hand side of (3.14) is well-defined 
since by (2.30), (2.38), and (3.11), 

[G p+1 (z,n) - y][G p+1 (z,n) + y] = G p+1 (z,n) 2 - R 2 p+2(z) 

= 4a(n) 2 F p {z,n)F+{z,n) 
v 

= 4a(n) 2 ]J[z- Mj (n)][z-^(n+l)] ^ 0, z£tt c . (3.15) 



Adding and subtracting g(z,n) on the right-hand side of (3.14) yields 

' - g{z,n)-\ — , z e u c , 



2dz \Gp +1 (z,n) + y 



y 



(3.16) 



where 



K{z,n) = -G p+ i(z,n) 



F;(z,n) + (F+)'(z, 



F p {z,n) F+(z,n) 
Next we prove that the mean value of K(z,-) equals zero. 
Lemma 3.4. Assume Hypothesis 3.2. Then 

(K(z,-))=0, zeIT c . 
Proof. Let z eU c . Using (2.28) we rewrite (3.17) as 



G; +1 (z,n) - F p (z,n). (3.17) 



(3.18) 



K(z,n) =-G p+1 (z,n) 



— In (G p+ i (z, n) + G p+1 (z,n)) 



+ d~z ln ( G P+i( z ' n ) + G f+i( z ' n )) 



Lr p+ i(z, n) + 



dz 



z — b(n) z — b + (n) 
G' p+1 {z,n) + (G; +1 )'(z,n) 

G p+1 (z,n) + G~ +1 (z,n) 

p+U '") ' 



+ 



(G+y(z,n) + G; +1 (z,n) 



Gp +1 (z,n) + G p+1 (z,n) 
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l{G p+1 (z,n) G p+1 {z,n) 



b{n) z-b+{n) 
(Gp +1 )'(z,n)G p+1 (z,n) - G' p+1 (z,n)G+ +1 (z,n) 
G+ +1 (z,n) + G p+1 (z,n) 
G p+1 {z,n)G p+1 {z,n) - (G p+1 )'{z,n)G p+ i{z,n) 



G p+ i(z,n) + G p+1 (z,ri) 

+ -n c . (3.19) 

2 \ z - b(n) z — b+(n) J 

Since K(z, •) is a sum of two difference expressions and G p+ i(z, •) and G* +1 (z, •) 
are bounded for fixed z G lie, one obtains 

(K{z,-))=0, zellc. (3.20) 

□ 

Using (3.16) and Lemma 3.4, one derives the following result that will subse- 
quently play a crucial role in this paper. 

Lemma 3.5. Assume Hypothesis 3.2 and let z,z G II. Then 



In 



Gp+i{z,-) -y 
G p+1 (z,-) + y 



' / dz'(g(z', 

J Zq 



)) 



In 



G p+ i(z , ■) - y 
G p+ i(z , -)+y 



(3.21) 



where the path connecting zq and z is assumed to lie in the cut plane II. Moreover, 
by taking limits to points on C in (3.21 ) ; the result (3.21) extends to either side of 
the cuts in the set C by continuity with respect to z. 

Proof. Let z, z G lie. Integrating equation (3.16) from z to z along a smooth 
path in lie yields 



hi 



G p+1 (z,n) - y 
G p +i(z,n) + y 



In 



Gp+i(z ,-) -y 
G p+ i(z a , -)+y 



! [ dz'g(z',n)+ 

J Zq 

,K(z',n) 



dz' 



By Lemma 3.3, K(z, •) is quasi-periodic. Consequently, also 

,K(z>,.) 



f 

J Zn 



dz' 



(3.22) 



(3.23) 



y 

is a family of uniformly almost periodic functions for z varying in compact subsets 
of lie as discussed in [22, Sect. 2.7]. By Lemma 3.4 one thus obtains 

dz' ^lil \ = 0. (3.24) 

y }/ 

Hence, taking mean values in (3.22) (taking into account (3.24)), proves (3.21) for 
z G lie- Since ft, I G No, are quasi-periodic by Lemma 3.3 (we recall that f = 1), 
(2.20) and (3.4) yield 



r dz'{g{z',-))=j2{f P - t ) r dz' 

J ZQ J Zq 



R 



2p+2 



(3.25) 



Thus, f* g dz' (g(z' , •)) has an analytic continuation with respect to z to all of II and 
consequently, (3.21) for z G Be extends by analytic continuation to z G B. By 
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continuity this extends to either side of the cuts in C. Interchanging the role of z 
and zo, analytic continuation with respect to zo then yields (3.21) for z, zo G II. □ 

Remark 3.6. For z G lie, the sequence m ( G P ^[z'!)+y ) is quasi-periodic and hence 



i(*v)-y 



)) is well-defined. But if one analytically continues 



ln ( %l\\\ Z z"n)+y ) with rcs P cct to z ' then (G P +i(z,n) - y) and (G p+1 (z,n) + y) may 
acquire zeros for some n G Z and hence ln( ^+^|~^ ) £ QP(Z). Nevertheless, 
as shown by the right-hand side of (3.21), admits an analytic 

continuation in z from lie to all of II, and from now on, ( m ( g P |^' z G II, 
always denotes that analytic continuation (cf. also (3.27)). 

Next, we will invoke the Baker- Akhiezer function ip(P,n,no) and analyze the 
expression (ln( g p +^'^~^ )) in more detail. 

Theorem 3.7. Assume Hypothesis 3.2, let P = (z,y) G n±, and n,no G Z. 

~ ~(3) ~(3) 

Moreover, select a homology basis {aj,bj}^ =1 on JC P such that B_ — f/ , with f/ 
the vector of b-periods of the normalized differential of the third kind, ui 
satisfies the constraint 



O) 



B = U 



(3) 



(cf. Appendix B). Then, 



Re < In 



Gp+ijz, -)-y 
G p +i(z, -) + y 



2Re 



'(3) 



Qo 



(3.26) 



(3.27) 



Proof. Using (2.41), (2.45) and (2.46) one obtains the following representation of 
the Baker- Akhiezer function tp(P, n, n ) for n > n , n, n n G Z, P G JC P , 



i>(P,n,n ) = [] ^ P ' m ) = 

m— no 

F p {z,n)\ 1/2 



y - G p+ i(z,m) F p (z,m+1) 



F p (z,n ) 
F P {z,n) 



n — 

-y- 

G p+1 (z,m) 



y-G p+ i(z,m) F p (z,m) 

1/2 

y 



- rn—no 
n-1 

n G p+ i(z,m) + y 

m=n ' 



1/2 



F p (z,n ) 



1/2 



x exp 



n-l 



1 " " 



In 



1 



x exp( -(n - n )( In 



G p+1 (z,m) - y 
G p+ i(z,m) + y 



G p+ i(z, -)-y 



- < In 



G p+ i(z, ■) -y 
G P +i(z,-) +y 



(3.28) 



G p+ i{z,-) + ?/, 
P = (z,y) G n±, z G U c , n,n G Z. 
A similar representation can be written for ip(P, n, uq) if n < no, n, no G Z, P G /C p . 



Since 



L \G p+ i(z,m)+y J 



' -t 71—1 r 
m—n 



In 



G p+ i(z,m) - ?y 
G p+ i(z,m) + y 



In 



has mean zero 
G p +i(z, •) -y 



G p+ i(^, •) +?y 



, = o(|n|), 
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z G Tic (3.29) 

by Theorem 3.1 (vi). In addition, the factor F p (z,n)/F p (z,no) in (3.28) is quasi- 
periodic and hence bounded on Z. 
On the other hand, (2.71) yields 



m n, no) = C(n, -o)^-^yy exp (^(n - n„) ^ ^ 

rP 



= e(P,n,n )exp \(n - n ) J ^pI + , Poo _J, (3-30) 

i , G/C p \({P 00± }U{A j (no)}? =1 ). 

Taking into account (2.67), (2.73), (2.81), (A.29), and the fact that by (2.60) no 
fij(n) can reach Poo ± as n varies in Z, one concludes that 

G(P,;n )e£°°(Z), P e IC p \{fi 3 (n )} P 3=1 . (3.31) 

A comparison of (3.28) and (3.30) then shows that the o(|n|)-term in (3.29) must 
actually be bounded on Z and hence the left-hand side of (3.29) is almost periodic 
(in fact, quasi-periodic). In addition, the term 



exp 



1 n — 1 

- y 

2 ^ 

m—no 



G p+1 (z,m) -y \ _L( Gp+i(z,-) -y 



G p+ i(z,m) + y J \ \G p+ i(z,-)+y 



z g n c , (3.32) 



is then almost periodic (in fact, quasi-periodic) by Theorem 3.1 (x). A further com- 
parison of (3.28) and (3.30) then yields (3.27) for z G He- Analytic continuation 
with respect to z then implies (3.27) for z e II. By continuity with respect to z, 
taking boundary values to either side of the cuts in the set C, this then extends to 
zeC (cf. (A.3), (A.4)) and hence proves (3.27) for P = (z,y) G /C P \{P 00± }. □ 



4. Spectra of Jacobi Operators with Quasi-Periodic 
Algebro-Geometric Coefficients 

In this section we establish the connection between the algebro-geometric for- 
malism of Section 2 and the spectral theoretic description of Jacobi operators H in 
f 2 (Z) with quasi-periodic algebro-geometric coefficients. In particular, we introduce 
the conditional stability set of H and prove our principal result, the characterization 
of the spectrum of H . Finally, we provide a qualitative description of the spectrum 
of H in terms of analytic spectral arcs. 

Suppose that a, b G £°°(Z) n QP(Z) satisfy the pth stationary Toda equation 
(2.18) on Z. The corresponding Jacobi operator H in £ 2 (Z) is then defined by 

H = aS + +a~S- + b, dom(ff ) = £ 2 (Z). (4.1) 

Thus, H is a bounded operator on £ 2 (Z) (it is self-adjoint if and only if a and b are 
real- valued) . 

Before we turn to the spectrum of H in the general non-self-adjoint case, we 
briefly mention the following result on the spectrum of H in the self-adjoint case 
with quasi-periodic (or almost periodic) real-valued coefficients a and b. We denote 
by cr(A), a e (A), and ad(A) the spectrum, essential spectrum, and discrete spectrum 
of a self-adjoint operator A in a complex Hilbert space, respectively. 
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Theorem 4.1 (See, e.g., [52] in the continuous context). Let a, b 6 QP(Z) be 
real-valued. Define the self-adjoint Jacobi operator H in £ 2 (Z) as in (4.1). Then, 

a(H) = a c (H) (4.2) 

C [-2sup(|Re(a(n))|) + inf (Re(6(n))) , 2 sup (|Re(a(n))| + supRe(6(n)))], 

nez ™e z nez nez 

<7d(-H") = 0- (4.3) 
Moreover, <j(H) contains no isolated points, that is, cr{H) is a perfect set. 

In the special periodic case where a, b are real- valued, the spectrum of H is 
purely absolutely continuous and a finite union of some compact intervals (see, 
e.g., [12], [13], [23], [53]-[56],). 

Next, we turn to the analysis of the generally non-self-adjoint operator H in 
(4.1). Assuming Hypothesis 3.2 we introduce the set £ C C by 

Re(U G ^-\- y )))=0Y (4-4) 



G p+ i{\-)+y / 

Below we will show that £ plays the role of the conditional stability set of H, famil- 
iar from the spectral theory of one-dimensional periodic differential and difference 
operators. 

Lemma 4.2. Assume Hypothesis 3.2. Then £ coincides with the conditional sta- 
bility set of H, that is, 

E = {A 6 C | there exists at least one bounded solution 

O^er(Z) of Hi, = Xtp}. (4.5) 

Proof. By (2.71) and (2.72), 

mn,n ) = C(n,n ) 9U{p - {no))) exp {(n - n ) J^^), (4.6) 

p = (z,y)en ± , 

is a solution of Hip = zip which is bounded on Z if and only if the exponential 
function in (4.6) is bounded on Z. By (3.27), the latter holds if and only if 

Re(U ^dzy )\)=0. (4.7) 



G p+1 (z,-) + y 



□ 



Remark 4.3. At first sight our a priori choice of cuts C for i? 2 p+2(-) 1 ^ 2 i as described 
in Appendix A, might seem unnatural as they completely ignore the actual spectrum 
of H. However, the spectrum of H is not known from the outset, and in the case 
of complex- valued potentials, spectral arcs of H may actually cross each other (cf. 
Theorem 4.7 (iv)) which renders them unsuitable for cuts of R,2 P +2(-) 1 ^ 2 ■ 

Before we state our first principal result on the spectrum of H, we find it con- 
venient to recall a number of basic definitions and well-known facts in connection 
with the spectral theory of non-self-adjoint operators (we refer to [18, Chs. I, III, 
IX], [27, Sects. 1, 21-23], [31, Sects. IV.5.6, V.3.2], and [49, p. 178-179] for more 
details). Let S be a densely defined closed operator in complex separable Hilbcrt 
space H. Denote by B(H) the Banach space of all bounded linear operators on 
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Ti and by ker(T) and ran(T) the kernel (null space) and range of a linear opera- 
tor T in TL. The resolvent set, p(S), spectrum, a(S), point spectrum (the set of 
eigenvalues), <J p (S), continuous spectrum, cr c (S), residual spectrum, cr r (5), field of 
regularity, ir(S), approximate point spectrum, <7 ap (5), two kinds of essential spec- 
tra, <r c (S), and a e (S), the numerical range of 5, 0(5), and the sets A (5) and A (5) 



arc defined as follows: 

p(S) ={z£C\(S-zI)- 1 €B(H)}, (4.8) 

a(S) = C\p(5), (4.9) 

<ip(5) = {A G C | kcr(5 - XI) + {0}}, (4.10) 
cr c (5) = {A G C | kcr(5 - XI) = {0} and ran(5 - XI) is dense in H 

but not equal to H}, (4-11) 
cr r (5) = {A G C | ker(5 - A/) = {0} and ran(5 - XI) is not dense in H}, (4.12) 
7r(5) = {z G C | there exists k z > s.t. ||(5 - z/)w|| w > 

for all u G dom(5)}, (4.13) 

a ap (S)=C\n(S), (4.14) 

A(5) = {z G C | dim(ker(5 - zl)) < oo and ran(5 - zl) is closed}, (4.15) 

<j e (5) = C\A(5), (4.16) 

A(5) = {z e C | dim(kcr(5 - zl)) < oo or dim(kcr(5* - z/)) < oo}, (4.17) 

a e {S) = C\A(5), (4.18) 

0(5) - {(/, 5/) e C | / e dom(S), ||/|| w - 1}, (4.19) 

respectively. One then has 

cr(5) = er p (5) U cr c (5) U cr r (5) (disjoint union) (4.20) 

= ( 7 p (5)Ua c (5)UCT r (5), (4.21) 

a c (S) C <7e(5)\(a p (5) U <j r (5)), (4.22) 

<j r (5) = <7 p (5*)*\a p (5), (4.23) 
c a p(5) = {A e C | there exists a sequence {f n }neN C dom(5) 

with ||/„|| w = 1, n e N, and lim ||(5 - A7)/„|| w = 0}, (4.24) 

n— >oo 

?c(5) C u e (5) C u ap (5) C ct(5) (all four sets are closed), (4-25) 

p(S) C tt(5) C A(5) C A(5) (all four sets are open), (4.26) 

a e (S) C 0(5), 0(5) is convex, (4.27) 

?c(5)=cr c (5) if 5 = 5*. (4.28) 

Here ct* in the context of (4.23) denotes the complex conjugate of the set a C C, 
that is, 

cr* = {A G C | A e cr}. (4.29) 



We note that there are several other versions of the concept of the essential spectrum 
in the non-self-adjoint context (cf. [18, Ch. IX]) but we will only use the two in 
(4.16) and in (4.18) in this paper. 

We start with the following elementary result. 
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Lemma 4.4. Let H be defined as in (4.1). Then, 

a c (H) = a c (H) C 6(F). (4.30) 

Proof. Since H and H * are second-order difference operators on Z, 

dim(kcr(ff - zl)) < 2, dim(kcr(iT - z/)) < 2. (4.31) 

Moreover, we note that S closed and densely defined and dim(ker(S'* — ~zl)) < oo 
implies that ran(S - zl) is closed (cf. [18, Theorem 1.3.2]). Equations (4.15)-(4.18) 
and (4.27) then prove (4.30). □ 

Theorem 4.5. Assume Hypothesis 3.2. Then the point spectrum and residual 
spectrum of H are empty and hence the spectrum of H is purely continuous, 

a p {H) = a r {H) = 0, (4.32) 
a(H) = <j c (H) = <j e (H) = a ap (H). (4.33) 

Proof. First we prove the absence of the point spectrum of H. Suppose z G II\(S U 
{pij(no)}P =1 ). Then ip(P, -,no) and ip(P*, -,«o) ar e linearly independent solutions 
of Hip = zip which are unbounded at +oo or — oo. This argument extends to all 
z G by multiplying tp(P,-,n ) and ip(P* , -,n ) with an appropriate function 
of z and no (independent of n). It also extends to either side of the cut C\S by 
continuity with respect to z. On the other hand, any solution ip(z, •) £ £ 2 (Z) of 
Hip = zip, z e C, is necessarily bounded (since any sequence in ^ 2 (Z) is bounded). 
Thus, 

{C\E}n(7 p (ff) = 0. (4.34) 

Hence, it remains to rule out eigenvalues located in S. We consider a fixed A 6 S 
and note that by (2.52), there exists at least one solution ^>i(A, ■) G ^°°(Z) of 
Hip — Aip. Actually, a comparison of (3.28) and (4.4) shows that one can choose 
V>i(A, •) such that |V>i(A, -)l G QP(Z) and hence Vi(A, •) £ £ 2 (Z). 

Next, suppose there exists a second solution ^(A, •) G £ 2 (Z) of Hip — Aip which is 
linearly independent of ipi(A, •). Then one concludes that the Wronskian of ^l^ ') 
and ip 2 (A,-) lies in £ 2 (Z), 

W(M\-),MK-)) el 2 {Z). (4.35) 

However, by hypothesis, W{ip\{A, •), ip 2 {A, ■)) = c{A) ^ is a nonzero constant. 
This contradiction proves that 

EnCTp(i7) = (4.36) 

and hence <y p (H) = 0. 

Next, we note that the same argument yields that H* also has no point spectrum, 



a p (H*) = 0. (4.37) 
Indeed, if a, b € £°°(Z) n QP(Z) satisfy the pth. stationary Toda equation (2.18) on 
Z, then a, b also satisfy one of the pth stationary Toda equation (2.18) associated 
with a hyperelliptic curve of genus p with {_B m }^Q replaced by {£ m }Si e ^ c - 
Since by general principles (cf. (4.29)), 

o-r(B) C a v {B*)* (4.38) 

for any densely defined closed linear operator B in some complex separable Hilbert 
space (see, e.g., [28, p. 71]), one obtains a r (H) = and hence (4.32). This proves 
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that the spectrum of H is purely continuous, a(H) = a c (H). The remaining 
equalities in (4.33) then follow from (4.22) and (4.25). □ 

The following result is a fundamental one: 

Theorem 4.6. Assume Hypothesis 3.2. Then the spectrum of H coincides with £ 
and hence equals the conditional stability set of H, 

= {A G C | there exists at least one bounded solution 

^ V G ^°°(^) of Hi) = AV>}. (4.40) 

In particular, 

{EmfZtl C o-(H), (4.41) 
and o-{H) contains no isolated points. 

Proof. First we will prove that 

a(H) C £ (4.42) 
by adapting a method due to Chisholm and Evcritt [10] (in the context of differential 
operators). For this purpose we temporarily choose z G n\(£ U {^j("o)}j = i) an d 
construct the resolvent of H as follows. Introducing the two branches ip±(P, n, n ) 
of the Baker- Akhiezer function if>(P, n, n ) by 

ip±(P,n,no) = i)(P,n,n ), P = (z, y) G IT±, n, no G Z, (4.43) 



we define 



)ip+(z,n,no) if "0+( z j ■> n o) G £ (no, oo), 
^+(z,n,n ) = ^ , , (4-44) 

ip-(z,n,no) if V-0v,™o) G £ z (n ,oo), 

V>-(z,n,n ) if -,n ) G ^ 2 (-oo,n ), 

ip + (z,n,n ) if t/> + (z, -,n ) G ^ 2 (-oo,n ), 
z G II\E, n, n G 



i>-(z,n,no) = { L ) ( ^ , )_ N ^ ( (4.45) 



and 

G(z, n, n') 



1 I ip-(z,n', n )tp+(z,n, no), n>n', 

W{ip-(z,n,n ),'4> + (z,n,n )) \i>-{z,n,no)i> + (z,n',n ), n<n', 

z G II\£, n,n G Z. (4.46) 

Due to the homogeneous nature of G, (4.46) extends to all z G II. Moreover, we 
extend (4.44)-(4.46) to either side of the cut C except at possible points in £ (i.e., 
to C\£) by continuity with respect to z, taking limits to C\£. Next, we introduce 
the operator R(z) in £ 2 (Z) defined by 

(R(z)f)(n) = ]T G(z,n,n')f(n'), f G £^(Z), z G II, (4.47) 

n'GZ 

where £g° (S) denotes the linear space of compactly supported (i.e., finite) complex- 
valued sequences, and extend it to z G C\£, as discussed in connection with 
G(-,n, n'). The explicit form of ip±(z, n, no), inferred from (3.30) by restricting 
P to n±, then yields the estimates 

\ijj±{z,n,n )\ < C±{z,n )e TK{z)n 7 z G II\£, n G Z (4.48) 
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for some constants C±(z,no) > 0, k(z) > 0, z G n\E. One can follow the second 
part of the proof of Theorem 5.3.2 in [20] line by line and prove that R(z), z G 
C\E, extends from £g°(Z) to a bounded linear operator defined on all of £ 2 (Z). A 
straightforward computation then proves 

(H - zI)R{z)f = f, f G £ 2 (Z), z G C\S (4.49) 

and hence also 

R{z)(H - zl)g = .9, g G £ 2 (Z), z G C\E. (4.50) 
Thus, R(z) = {H- zl)' 1 , z G C\E, and hence (4.42) holds. 
Next we will prove that 

a{H) D E. (4.51) 

We will adapt a strategy of proof applied by Eastham in the continuous case of (real- 
valued) periodic potentials [19] (reproduced in the proof of Theorem 5.3.2 of [20]) 
to the (complex-valued) quasi-periodic discrete case at hand. Suppose A G E. By 
the characterization (4.5) of E, there exists a bounded solution •) of Hip = 
A comparison with the Baker-Akhiezer function (3.30) then shows that one can 
assume, without loss of generality, that 

MV)|eQP(Z). (4.52) 

By Theorem 3.1 (i), one obtains 

V>(A,-) Gf°(Z). (4.53) 
Next, we pick O G N and consider g(n), n = 0, 1, . . . , fl, satisfying 

5(0) = 0, g (n) = i, 
0<g(n)<l, n = 
Moreover, we introduce the sequence {/ifc}fceN G ^ 2 (Z) by 

1, |n|<(fc-l)ft, 
h k (n) = I g(kQ - \n\), (k - < \n\ < fcft, (4.55) 
0, \n\ > kfl 

and the sequence {fkWjkeN G £ 2 (Z) by 

/fe(A, n) = d k (X)i>(\, n)h k (n), neZ, 4(A) > 0, k G N. (4.56) 
Here d k (A) is determined by the normalization requirement 

||/ fc (A)|| 2 = l, keN. (4.57) 

Of course, 

f k (X,-)e£ 2 (Z), fcGN, (4.58) 

since /fe(A, •) is finitely supported. Next, we note that as a consequence of Theorem 
3.1 (viii), 



(4.54) 



N 

5>(A,n)| 2 = (2A + l)(|^(A,.)| 2 )+o(A) (4.59) 

-JV 

with 

2\ 



<l^(A,-)| 2 )>0. (4.60) 
Thus, one computes 

l = ll/ fe (A)||^ = 4(A) 2 ^|^(A,n)| 2 ^(n) 2 
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= 4(A) 2 \^{\n)\ 2 h k {nf>d k {\) 2 Y, l^( A ' n )! 2 

\n\<ki1 |n|<(fc-l)fi 

> 4(A) 2 [<|V(A, -)| 2 )(fc - 1)0 + o(k)] . (4.61) 

Consequently, 

d k (X) = Oik- 1 / 2 ). (4.62) 

k— >oo 

Next, one computes 

(if - A7)/ fe (A, n) -4(A) [a(n)V(A, n) [h k (n + 1) - /i fe (n)] 

+ a(n - 1)V(A, n - 1) [/i fe (n - 1) - /*(»)]] (4.63) 

and hence 

||( J ff-A/)M| 2 <24(A)||a|| 0O ||^(A)(/ l +-4)|| 2 , fceN. (4.64) 
Using (4.53) and (4.55) one estimates 

M\)[h+ -h k ]\\ 2 2 = Yl m\n)\ 2 \h k (n + l)-h k (n)\ 2 

(fc-l)0<|n|<fcO 

< 211^)^(0 + 1). (4.65) 
Thus, combining (4.62) and (4.64)-(4.65) one infers 

lim ||(tf-A/)/ fe || 2 = (4.66) 

n^oo 

and hence A G a &p (H) = a(H) by (4.24) and (4.33). 

Relation (4.41) follows from (4.5) and the fact that by (2.52) there exists a 
solution ip((E m , 0), •, n ) G £°°(Z) of Hip = E m ip for all m = 0, . . . , 2p + 1. 

Finally, u(H) contains no isolated points since those would necessarily be es- 
sential singularities of the resolvent of H, as H has no eigenvalues by (4.32) (cf. 
[31, Sect. III. 6. 5]). An analysis of the Green's function of H reveals at most an 
algebraic singularity at the points {Em} 2 ^^ and hence excludes the possibility of 
an essential singularity of (H — zl)~ x . □ 

In the special self-adjoint case where a, b are real-valued, the result (4.39) is 
equivalent to the vanishing of the Lyapunov exponent of H which characterizes the 
(purely absolutely continous) spectrum of H as discussed by Carmona and Lacroix 
[9, Chs. IV, VII] (cf. also [8], [14], [29], [32]). 

The explicit formula for £ in (4.4) permits a qualitative description of the spec- 
trum of H as follows. We recall (3.16) and (3.25) and write 

IdA /g^Q-yNV n^*-*'* z G II, (4.67) 

for some constants 

feiCC. (4.68) 

As in similar situations before, (4.67) extends to either side of the cuts in C by 
continuity with respect to z. 
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Theorem 4.7. Assume Hypothesis 3.2. Then the spectrum cr(H) of H has the 

following properties: 

(i) a(H) C C is bounded, 



a(H) C {z G C | Re(z) G [M 1; M 2 ], Im(z) £ [M 3 , M 4 }}, (4.69) 



where 



M x = -2sup[|Re(a(n))|] + inf [Re(6(n))], 
M 2 = 2sup[|Re(a(n))|] +sup[Re(6(n))], 

Af 3 = -2sup[|Im(a(n))|] + inf [Im(6(n))], 

nez ™ez 

M 4 = 2sup[|Im(o(n))|] + sup[Im(6(n))]. 

nGZ neZ 



(4.70) 



(m) cr(H) consists of finitely many simple analytic arcs (cf. Remark 4-8)- These 
analytic arcs may only end at the points Ai, . . . , X p , Eq, . . . , E2 P +i- 
(Hi) Each E m , m = 0, . . . , 2p+l, is met by at least one of these arcs. More precisely, 
a particular E mo is hit by precisely 2N$ + 1 analytic arcs, where Nq e {0,...,p} 
denotes the number of Xj that coincide with E mo . Adjacent arcs meet at an angle 
2tt/(2Nq + 1) at E ma . (Thus, generically, Nq = and precisely one arc hits E mo .) 
(iv) Crossings of spectral arcs are permitted. This phenomenon takes place precisely 
when for a particular j G {1, . . . ,p}, Xj G cr(H) such that 

Re I ( ln{ Gp+1 ^ j0 '' ) - y \ ) ) =0 



G p+1 (X j0 ,-)+yJ I ) (4.71) 
for some j G {1, . . . ,p} with X jo £ {E m }^Q. 

In this case 2Mq-\-2 analytic arcs are converging toward Xj , where Mq G {1, . . . , p} 
denotes the number of Xj that coincide with Xj . Adjacent arcs meet at an angle 
7r/(M + 1) at Xj . (Thus, if crossings occur, generically, M = 1 and two arcs 
cross at a right angle.) 

(v) The resolvent set C\a(H) of H is path- connected. 

Proof. Item (i) follows from (4.30) and (4.33) upon noticing that 

oo 

(f,Hf) = 2 a(k)Re[f(k + l)J(k)} + (f,Re(b)f)+i(f,lm(b)f), f G l 2 (Z). 

k— — oo 

(4.72) 

To prove (ii) we first introduce the meromorphic differential of the third kind 

y R2 P +2(z) i/2 

P=(z,y)elC p \{Po 0± } (4.73) 

(cf. (4.68)). Then, by Lemma 3.5, 

\ \G p+1 (z,-)+yJ / J Qo \ \G p+1 (z ,-)+yJ /' 
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for some fixed Qo = (zo,yo) £ ICp\{Poo±}, is holomorphic on IC p \{Poc ± }- By 
(4.67), (4.68), the characterization (4.39) of the spectrum, 



a(H) = Ae 



and the fact that Re(( In ( g p **[*'')+^ ))) is a harmonic function on the cut plane IT, 
the spectrum a(H) of H consists of analytic arcs which may only end at the points 
Ai, . . . , A p , Eq, . . . , E2 P +\- (Since cr(H) is independent of the chosen set of cuts, if 
a spectral arc crosses or runs along a part of one of the cuts in C, one can slightly 
deform the original set of cuts to extend an analytic arc along or across such an 
original cut.) 

To prove (Hi) one first recalls that by Theorem 4.6 the spectrum of H contains 
no isolated points. On the other hand, since {i? m }m=o c a (H) by (4.41), one 
concludes that at least one spectral arc meets each E m , m = 0, . . . , 2p+l. Choosing 
Qo = (-Emo)O) in (4.74) one obtains 

m ( G jr^-\~ yS \ )=*f •)) + ( m f G ^ (E " 



G p+ i(z,-)+yJ / J Emo \ \G p+ i(E, 

of dz' fej^j / ( Gp+1 (E mo ,-)-y 
J^o m^'-E m )) 1/2 \ \G p+ i(E mo ,-) + y 

f dz' (z' - E mo ) N "-^ [C + 0(z' - E ma )] 

o JE mn 

G p +i(E mo , •) — y 





In 



G p+ i(E mo ,-) + y 

(z E mo ) N °+(W [c + E ^ }] + / / G +l (E ■ ) jf (4 __ (i) 



z^E mo JV + (l/2) 1 v man \ \G p+1 (E mo ,-)+y 
for some C = \C\e^° G C\{0}. Using 

Re // ln ( G p+1 (E m ,-)-y \ )]=() m = () , p -| (4 _ 77) 



G p+ \(E m , ■) + y 

as a consequence of (4.41), Re(( ln = and z = E mo + pe** imply 

= cos[(7V + (1/2)) V + ^o]p Ar ° +(1/2) [|C| + 0(p)]. (4.78) 
pio 

This proves the assertions made in item (Hi). 

In order to prove (iv) it suffices to refer to (4.67) and observe that locally 

ls( ln ( Gptifc-)+y )) behaves like °o( z - K) Mo for some C £ C\{0} in a suf- 
ficiently small neighborhood of Aj . 

Finally we will show that all arcs are simple (i.e., do not self- intersect each other). 
Assume that the spectrum of H contains a simple closed loop 7, 7 C cr(H). Then 
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where the closed simple curve r C IC P denotes an appropriate lift of 7 to IC P , yields 
the contradiction 

Re (7 In f ^ P+1 ! Z S^?''?~ y j^? > ) \] = for all P in the interior of T (4.80) 

by Corollary 8.2.5 in [3]. Therefore, since there are no closed loops in <j(H) and 
no analytic arc tends to infinity, the resolvent set of H is connected and hence 
path-connected, proving (v). □ 

Remark 4.8. Here a C C is called an arc if there exists a parameterization 7 G 
C([0, 1]) such that a = {"f(t) \ t G [0, 1]}. The arc a is called simple if there exists a 
parameterization 7 such that 7: [0, 1] — > C is injective. 

Appendix A. Hyperelliptic Curves and their Theta Functions 

We provide a brief summary of some of the fundamental notations needed from 
the theory of hyperelliptic Riemann surfaces. More details can be found in some 
of the standard textbooks [21] and [44] as well as in monographs dedicated to 
integrable systems such as [4, Ch. 2], [24, App. A, B]. 

Fix peN. We intend to describe the hyperelliptic Riemann surface K p of genus 
p of the Toda-type curve (2.36), associated with the polynomial 

F P (z,y) = y 2 - i? 2p +2(z) = 0, 

2p+l (A 1) 

R2 P+ 2(Z) = 1] ( Z E ™\ {EmfZtl C C. 
m— 

To simplify the discussion we will assume that the affine part of K, p is nonsingular, 
that is, we assume that 

E m ^ E m t for m^m', m, m' = 0, . . . , 2p + 1 (A. 2) 

throughout this appendix. Next we introduce an appropriate set of (nonintersect- 
ing) cuts Cj joining E m ^ and E m i^, j = 1, . . . ,p + 1, and denote 

P+i 

C=[jCj, C,nC fc = 0, j^k. (A.3) 
j=i 

Define the cut plane 

II = C\C, (A.4) 
and introduce the holomorphic function 

- 2p+ l N 1/2 

(A.5) 



('FT' \ 
n ( z e ^ ) 



on II with an appropriate choice of the square root branch in (A.5). Next we define 
the set 

M p = {(z,aR 2p+2 (z) 1 / 2 ) \zeC, <J G {1, 1}}U {^00+ > Poo- } (A.6) 

by extending i?2p+2(') 1 ^ 2 to C. The hyperelliptic curve JC P is then the set M. p with 
its natural complex structure obtained upon gluing the two sheets of A4 p crosswise 
along the cuts. Moreover, we introduce the set of branch points 

B(K p ) = {(E m ,Q)}^tl (A.7) 
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Points P G /C p \{Poo ± } are denoted by 

P=(z,aR 2p+2 (z) 1 / 2 ) = (z,y), (A.8) 

where y{P) denotes the meromorphic function on K p satisfying T p {z,y) = y 2 — 
R2 P +2(z) = and 

V(P) ( :,t(iJ(E^)( + 0(C 2 )) C^ 1 as P - P oc± , C = l/z. (A.9) 

In addition, we introduce the holomorphic sheet exchange map (involution) 
*:/C p ^/C p , P=(z,y)~P* = (z,-y), P oc± ^P^ ± =P OOT (A.10) 
and the two meromorphic projection maps 

7r: K,p — > C U {oo}, P = (z, y) ^ z, Poo ± >— ► °o (A. 11) 

and 

y: Kp -^CU{oo}, P = (z, y) h-> y, P QO± ^ oo. (A.12) 
Thus the map 7r has a pole of order 1 at P QO± and y has a pole of order p + 1 at 
P oc± . Moreover, 

ff(P*) = 7f(P), y(P*) = -y(P), Pg/C p . (A.13) 

As a result, K, p is a two-sheeted branched covering of the Riemann sphere CP 1 
(=CU {oo}) branched at the 2p + 4 points {(E m , O)}^^ 1 , Poo ± ■ fC p is compact 
since 7f is open and CP 1 is compact. Therefore, the compact hyperelliptic Riemann 
surface resulting in this manner has topological genus p. 
Next we introduce the upper and lower sheets H± by 

n± - {(z, ±R 2p+2 (z) 1 i 2 ) e M P | z € n} (A.u) 

and the associated charts 

C±:il ± ^n, p^ z . (A.15) 

Let {dj, bj}^ =1 be a homology basis for JC P with intersection matrix of the cycles 
satisfying 

a,job k =Sj,k, a,joa k = 0, bjob k = 0, j,k = l,...,p. (A. 16) 

Associated with the homology basis {cij,bj}^ =1 we also recall the canonical dis- 
section of JC P along its cycles yielding the simply connected interior K, p of the 
fundamental polygon dK. p given by 

dlC p = a 1 b 1 ai 1 bi 1 a 2 b2<h 1 f>2 1 ' ' ' a^b' 1 . (A.17) 

Let M(IC P ) and A4 1 (JC P ) denote the set of meromorphic functions (0- forms) and 
meromorphic differentials (1-forms) on JC P , respectively. The residue of a meromor- 
phic differential v G A4 1 (JC P ) at a point Q G JC P is defined by 

reS ^ = i/ 7 /' (A ' 18) 

where 7q is a counterclockwise oriented smooth simple closed contour encircling Q 
but no other pole of v. Holomorphic differentials are also called Abelian differentials 
of the first kind. Abelian differentials of the second kind lo^ G A4 1 (/C p ) are 
characterized by the property that all their residues vanish. Any meromorphic 
differential lu^ on K, p not of the first or second kind is said to be of the third kind. 
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A differential of the third kind u/ 3 ) G M 1 (JC P ) is usually normalized by vanishing 
of its a-pcriods, that is, 



/ 

J a-i 



u>W=0, j = l,...,p. (A.19) 

l aj 
(3") 

A normal differential uj p p , associated with two distinct points P\, P2 G 1C P , by 
definition, has simple poles at Pi and Pi with residues +1 at Pi and —1 at P 2 and 

(3) 

vanishing a-periods. If u P q is a normal differential of the third kind associated 
with P, Q G K p , holomorphic on K, p \{P, Q}, then 

/rQ 
4% = 27ri J p lu 3 , j = l,..., p. (A.20) 

We shall always assume (without loss of generality) that all poles of a/ 3 ) on K p lie 
on IC P (i.e., not on dlC p ). 

Using local charts one infers that dz/y is a holomorphic differential on K p with 
zeros of order p — 1 at Poo ± and hence 

z^ 1 dz , , . 

Vj = » j = l,...,p, A.21 

y 

form a basis for the space of holomorphic differentials on K, p . Introducing the 
invcrtiblc matrix C in C p 



C =( C ^) jk=1 ... p , Cj, k = / Vj , (A.22) 



a,, 



c(k) = ( Cl (k),...,c p (k)), c j (k) = (C- 1 ) , j,k = l,...,p, (A.23) 



the normalized differentials ujj for j = 1, . . . ,p, 

P r 
ujj = c 3 (t)Vi, / = 5 j>k , j, k = 1, . . . ,p, 

form a canonical basis for the space of holomorphic differentials on JC p . 
In the chart {Up oo± , Cpoo ± ) induced by 1/n near P oc± one infers, 



(A.24) 



c(j)C p "^C 

,=i(nro(i-c^)) 
/ i 2p+1 \ 

= ± c(p) + ([-c(p) V P m + c(p - 1)] + 0(C 2 ) W as P -> Poc ± , 

C = 1/z. 

The matrix r = ( T j,f)^_ 1 in C pxp of 6-periods defined by 

Tj,e= [ u t , j,£ = l,...,p (A.26) 



Jb 3 

satisfies 

Im(r) > and tjj = n tj , j, I = 1, . . . , p. (A. 27) 

Associated with the matrix r one introduces the period lattice 

L p = {z£C p I z = m + nr, m,nEZ p } (A.28) 
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and the Riemann theta function associated with K p and the given homology basis 
{aj, 6j}j=i,...,p, 

= cxp ( 27 ™fe z) + Ki{n, nr)) , zeC p , (A.29) 

where (u,v) = m t = X^=i% v j denotes the scalar product in C p . It has the 
following fundamental properties 

0(zi, . . .,Zj-!, -Zj,z j+1 , ...,z n )= 9(z), (A.30) 

6{z + JR + H T ) = CX P ( — 2ni(n, z) — ni(n, nr))6(z_) 7 m,n^Z p . (A. 31) 

Next we briefly describe some consequences of a change of homology basis. Let 

{a 1 ,...,a p ,b 1 ,...,b p } (A.32) 

be a canonical homology basis on K p with intersection matrix satisfying (A. 16) and 
let 

{o , 1 ,...,o;,6i,...,&;} (A.33) 
be a homology basis on JC p related to (A.32) by 

f -^=xK), (AM) 



b' 



where 



= (ai, . . . ,a p ) T , 6 T = 0i, . . . ,b p ) T , 
= K,...,a p ) T , 6' T = (b[,..., b' p ) T , (A.35) 

x ={c d) ( a - 36 ) 

with A, B, C, and Z? being p x p matrices with integer entries. Then (A.33) is also 
a canonical homology basis on IC P with intersection matrix satisfying (A. 16) if and 
only if 

AeS P (p,Z), (A.37) 

where 

W) = {x=(£ B D ) a(_° /p £)*T=(_° Jp ; p )' dct(XH; 

(A.38) 

denotes the symplectic modular group (here A, B, C,D in X are again pxp matrices 
with integer entries). If {u)j} P =1 and {oj'j} P =1 are the normalized bases of holomor- 
phic differentials corresponding to the canonical homology bases (A.32) and (A.33), 
with r and r' the associated b and b'-pcriods of u = ui, . . . , u> p and u/ = u>[, . . . , ui' p , 
respectively, then one computes 

u/ = u(A + Bt)- 1 , t' = (C + Dt)(A + Bt)- 1 . (A.39) 

Fixing a base point Qo £ ^p\{Poo ± }, one denotes by J(IC P ) = C p /L p the Jacobi 
variety of JC P , and defines the Abel map Aq q by 

A Qo - IC n ^ J{K. P ), A Qa (P)=( / uji,..., ujp) (mod L p ), P e K p . 

\JQo JQo / 

(A.40) 
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Next, consider the vector of 6-pcriods of ujp^ p ^ /(2ni), the normalized 
differential of the third kind, holomorphic on /C P \{P 00± }, 

Using (A. 20) one then computes 

U 3) = A Poo _ (P 00+ ) = 2A Qa (P OQ+ ), (A.42) 

where Qo is chosen to be a branch point of K p , Qo E B(K, p ), in the last part of 
(A.42). 

Similarly, one introduces 

a QQ : Div(/C p ) - J(/C p ), P ~ a Qo (V) = £ V(P)A Qa (P), (A.43) 

where Div(A^ p ) denotes the set of divisors on K, p . Here a map P: fC p — > Z is called 
a divisor on /C p if P(P) 7^ for only finitely many P E JC P . (In the main body of 
this paper we will choose Qo to be one of the branch points, i.e., Qo G B(IC P ), and 
we will always choose the same path of integration from Qo to P in all Abelian 
integrals.) For subsequent use in Remark A. 4 we also introduce 



'Qo JQo 



Aq ■■ 1C P - U> , (A.44) 

P~A Qo (P) = {A QoA (P),...,A Qo!P (P)) = ( I wi, 

\JQo 

and 

S Qo : Div(/C p ) - C p , P ^ a Qo (P) = £ P(P)l Qo (P). (A.45) 

In connection with divisors on £ p we will employ the following (additive) nota- 
tion, 

Vq q = P Qo + Vq , Vq = V Qi +---+ V Qm , (A.46) 
Q = {Qi, . . . , Q m } G Sym" 1 /C p , Qo G /C p , m G N, 
where for any Q G /C p , 

and Sym m /C p denotes the mth symmetric product of IC p . In particular, Sym™ 1 K p 
can be identified with the set of nonnegative divisors < P G Div(/C p ) of degree 
m G N. 

For / G A4(/Cp)\{0}, w G A / t 1 (/C p )\{0} the divisors of / and lu are denoted 
by (/) and (w), respectively. Two divisors P, £ G Div(/C p ) are called equivalent, 
denoted by P ~ £, if and only if V-£ = (/) for some / G M(IC p )\{0}. The divisor 
class [P] of P is then given by [P] = {£ E Div(/C p ) | £ ~ P}. We recall that 

deg((/)) = 0, deg((w)) - 2(p-l), / G M(/C p )\{0}, w G M 1 (/C p )\{0}, (A.48) 

where the degree deg(P) of P is given by deg(P) = J2peK p ^(P)- ^ ^ s customary 
to call (/) (respectively, (w)) a principal (respectively, canonical) divisor. 
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Introducing the complex linear spaces 

C{V) = {/ e M{K P ) | / = or (/) > V}, r{V) = dim c C(V), (A.49) 

C\V) = {uj G M^/Cp) | w = or (w) > ©}, i(X>) - dim c &(V) (A.50) 

(with the index of specialty of T>), one infers that deg(X>), r(D), and i(T>) only 
depend on the divisor class [V] of V. Moreover, we recall the following fundamental 
facts. 

Theorem A.l. Let V e Div(/C p ), w £ A1 1 (/C P )\{0}. T/ien, 

t(X>) = r(D - H), peN . (A.51) 

The Riemann-Roch theorem reads 

r{-V) =deg(D)+i(D)-p+l, n e N . (A.52) 

-By Abel's theorem, V e Div(/C p ), peN, is principal if and only if 

deg(P) = and a Qo (£>) = 0. (A.53) 

Finally, assume p e N. TTien aQ : Div(/C p ) — ► J(fC p ) is surjective (Jacobi's 
inversion theorem). 

Theorem A.2. Let V Q e Sym p £ p , Q = {Qi, . . . , Q p }. T/ien, 

1 < %{Vg) = s (A.54) 

if and only if there are s pairs of the type {P, P*} C . . . , Q p } (this includes, 
of course, branch points for which P = P*). One has s < p/2. 

Next, we denote by Eq o — (£q 0j1 , . . . , Sq p ) the vector of Riemann constants, 

1 P f f P 

2q„, 3 = x(l + ^)-]r / u e {P) / Wj) i = l,...,p. (A.55) 

Theorem A. 3. Let Q = {Qi, . . . , Q p } £ Sym p K p and assume Vq to be nonspecial, 
that is, i{T>Q) = 0. Then, 

0(z Qo - A Qo (P) + <*«„(%)) = if and onl V ifPe{Qi,...,Qp}. (A.56) 

Remark A. 4. In Section 2 we dealt with theta function expressions of the type 

_ 0£ Qa A Qo (P) +a q „(P 1 )) / /-" A 
^ ( ) "^ 0o -A Qo( P) + a Qo( P 2 )) eXP v J' PG/C - (A ' 57) 

where 2?j e Sym p /C p , j = 1,2, are nonspecial positive divisors of degree p, c e C is 
a constant, and fi'- 3 ) is a normalized differential of the third kind with a prescribed 
singularity at P oc± . Even though we agree to always choose identical paths of 
integration from P to P in all Abelian integrals (A.57), this is not sufficient to 
render tp single- valued on /C p . To achieve single- valuedness one needs to replace /C p 
by its simply connected canonical dissection IC P and then replace Aq and Oq in 
(A.57) with Aq o and Sq o as introduced in (A. 44) and (A. 45). In particular, one 
regards dj,bj, j — 1, . . . ,p, as curves (being a part of d!C p , cf. (A. 17)) and not 
as homology classes. Similarly, one then replaces Sq o by Eq o (replacing Aq o by 
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Aq m (A. 55), etc.). Moreover, in connection with ip, one introduces the vector of 
6-periods U {3) of Q (3 *> by 

= &*>,... ,U<»), uf^^-J^V, j = l,...,p, (A.58) 

and then renders ip single-valued on JC P by requiring 

a Qo (V 1 )-a Qo (V 2 ) = cU^ (A.59) 

(as opposed to merely aQ g (Di) — aQ (T> 2 ) = cL^ 3 ** (mod L p )). Actually, by (A. 31), 

S Qo (Pi) -a Qo (P 2 ) - cU^ e TP, (A.60) 

suffices to guarantee single- valuedness of tp on K p . Without the replacement of 
Aq and Oq o by Aq o and Sq in (A. 57) and without the assumption (A.59) (or 
(A.60)), V is a multiplicative (multi-valued) function on fC p , and then most effec- 
tively discussed by introducing the notion of characters on K p (cf. [21, Sect. III. 9]). 
For simplicity, we decided to avoid the latter possibility and throughout this paper 
will always tacitly assume (A.59) or (A.60). 

Appendix B. Restrictions onB = C/^ 3) 

The purpose of this appendix is to prove the result (2.81), B = ' € W, for 
some choice of homology basis {aj, bj} P =1 on IC P as recorded in Remark 2.7. 

To this end we first recall a few notions in connection with periodic meromorphic 
functions of p complex variables. 

Definition B.l. Let p G N and F: C p — > C U {oo} be meromorphic (i.e., a ratio 
of two entire functions of p complex variables) . Then, 

(i) k> = ( w i) • • ■ 7 w p) £ C p \{0} is called a period of F if 

F(z + w)=F(z) (B.l) 

for all z£C p for which F is analytic. The set of all periods of F is denoted by Vf- 

(ii) F is called degenerate if it depends on less than p complex variables; otherwise, 
F is called nondegenerate. 

Theorem B.2. Let p E N, F: C p — > C U {oo} be meromorphic, and Vf be the set 
of all periods of F . Then either 

(i) Vf has a finite limit point, 
or 

(ii) Vf has no finite limit point. 

In case (i), Vf contains infinitesimal periods (i.e., sequences of nonzero periods 
converging to zero). In addition, in case (i) each period is a limit point of periods 
and hence Vf is a perfect set. 

Moreover, F is degenerate if and only if F admits infinitesimal periods. In partic- 
ular, for nondegenerate functions F only alternative (ii) applies. 

Next, let uj^ g C p \{0}, q = 1, . . . , r for some reR Then lu_ 1i . . . , w r are called 
linearly independent over Z (resp. R) if 

v\w_\ + • • • + f r w r = 0, v q e Z (resp., v q 6R), q = 1, . . . , r, 

implies v\ = ■ ■ ■ = v r = 0. (B-2) 
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Clearly, the maximal number of vectors in C p linearly independent over R equals 
2p. 

Theorem B.3. Let peN. 

(i) If F : C p — > C U {00} is a nondegenerate meromorphic function with periods 
Lo q G C p \{0}, q = 1, . . . ,r, r G N, linearly independent over Z 7 then Wj, . . . ,w r are 
a?so linearly independent over K. Ira particular, r < 2p. 

(ii) A nondegenerate entire function F : C p — > C cannot have more than p periods 
linearly independent over Z (orR). 

For p — 1, exp(z), sin(^) are examples of entire functions with precisely one 
period. Any non-constant doubly periodic meromorphic function of one complex 
variable is elliptic (and hence indeed has poles). 

Definition B.4. Let p,r E N. A system of periods uj q G C p \{0}, q = 1, . . . , r, of a 
nondegenerate meromorphic function F : C p — > CU{oo}, linearly independent over 
Z, is called fundamental or a 6asis of periods for F if every period w of F is of the 
form 

lu = miWj + • • • + 77j r w r for some m q G Z, o = 1, . . . , r. (B-3) 

The representation of w in (B.3) is unique since by hypothesis w 1; . . . , w r are 
linearly independent over Z. In addition, Vf is countable in this case. (This rules 
out case (?) in Theorem B.2 since a perfect set is uncountable. Hence, one does not 
have to assume that F is nondegenerate in Definition B.4.) 

This material is standard and can be found, for instance, in [39, Ch. 2]. 

Next, returning to the Riemann theta function 0(J in (A. 29), we introduce the 
vectors {e^ =11 { Zj } p j=1 C C p \{0} by 

e j = (0,...,0,^l > ,0,...,0), T 3 =e 3 T, j = l,...,p. (B.4) 

3 

Then 

{e 3 }U (B-5) 
is a basis of periods for the entire (nondegenerate) function 0(J : C p — > C. Moreover, 
fixing k G {1, . . . then 

!^-Z,^ ; (B.6) 

is a basis of periods for the meromorphic function d Zk In (g^pyy) : C p — > C U {00}, 
V G C p (cf. (A.31) and [21, p. 91]). 

Next, let A G C p , D = (£>i, . . . , D p ) G M p , G M\{0}, j = and 
consider 



/ fe :R^C, / fc (n) = a afc In 



9^ In 



0(A - z) 
0(C - z) 

0{A- zdmg{D))\ 



z—Dn 



(B.7) 

z=(n,...,n) 



^(C-zdiag^)) J 
Here diag(Z2) denotes the diagonal matrix 

diagGD) = , r (B.8) 
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Then the quasi-periods DJ 1 , j = 1, • • • ,P , of fk arc in a one-to-one correspondence 
with the periods of 

of the special type 

^•(diag(D)) _1 = (0,...,0,£>7 1 ,0,...,0). (B.10) 

Moreover, 

/fc(n)=F fc (2)|, =(n ,... in) , neZ. (B.ll) 

Theorem B.5. Suppose a and b in (2.75) to be quasi-periodic. Then there exists 

~ ~(3) ~(3) 

a homology basis {a.j,bj^ v - =1 on JC P such that the vector B_ = U_ with U_ the 
vector of b-periods of the corresponding normalized differential of the third kind, 
Wp2 + p 00 _; satisfies the constraint 

~ ~C31 

B = Uo eR p . (B.12) 

Proof. By (A. 41), the vector of 6-periods associated with a given homology 
basis {cij, 6j}j = i on JC P and the normalized differential of the third kind, wj, 3 ' p _ , 
is continuous with respect to i?o, • • • , ^2p+i- Hence, we may assume in the following 
that 

Bj ^ 0, .7 = 1,... ,p, B = (Bi, . . . ,B P ) (B.13) 
by slightly altering .Eg, ■ • ■ , .E^p+i, if necessary. Using (2.76), we may write 



w=0 

(B.14) 



z=Bn 



B = U^ ) . (B.15) 



where by (2.78), 

Introducing the meromorphic (nondegenerate) function V: C p — > C U {oo} by 

vU)^-t^-( ^I^i ). 

one observes that 

&(n)=VG0U=(„,...,„). (B.17) 
In addition, V has a basis of periods 

{ej ( diag(B)) ^ , Zj ( diag(fl)) _1 }* =i (B.18) 

by (B.6), where 

^.(diagGB))- 1 = (0,...,0,B7 1 ,0,...,0), j = l,...,p, (B.19) 

3 

T^diag^))" 1 = (T j , 1 B^ 1 ,...,T j!P B- 1 ), j = l,...,p. (B.20) 
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By hypothesis, b in (B.14) is quasi-periodic and hence has p real (scalar) quasi- 
periods. The latter are not necessarily linearly independent over Q from the out- 
set, but by slightly changing the locations of branchpoints {£ m }^ into, say, 
{E m }m=o , one can assume they are. In particular, since the period vectors in (B.18) 
are linearly independent and the (scalar) quasi-periods of b are in a one-one corre- 
spondence with vector periods of V of the special form (B.19) (cf. (B.9), (B.10)), 

~ ~(3) 

there exists a homology basis {ctjibj}^ =1 on fC p such that the vector B_ = U_ , 
corresponding to the normalized differential of the third kind, ui p p and this 

~ 3 

particular homology basis, is real-valued. By continuity of U_ with respect to 
Eq, . . . , E 2p +i, this proves (B.12). □ 

Acknowledgments. We are grateful to Leonid Golinskii for pointing out ref- 
erences [45] and [46] to us. 
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